Astronomy & Astrophysics manuscript no. mir'survey2'astroph 


© ESO 2010 


March 1, 2010 





The AKARI/IRC Mid- Infrared All-Sky Survey 

Daisuke Ishihara^'^, Takashi Onaka^, Hirokazu Kataza'^, Alberto Salama^, Carlos Alfageme^ *, Angelo 
Cassatella^'^'^, Nick Cox^ **, Pedro Garcia-Lario'^ , Craig Stephenson^ Martin Cohen''', Naofumi 
Fujishiro^'^ ^, Hideaki Fujiwara^, Sunao Hasegawa^, Yoshifusa Ita^, Woojung Kim^'^ ^, Hideo Matsuhara'^, 
Hiroshi Murakami^, Thomas G. Miiller^^, Takao Nakagawa^, Youichi Ohyama^^, Shinki Oyabu^, Jeonghyun 
Pyo^^, Itsuki Sakon^, Hiroshi Shibai^^, Satoshi Takita^, Toshihiko Tanabe^^, Kazunori Uemizu^, Munetaka 
Ueno^, Fumihiko Usui^, Takehiko Wada^, Hidenori Watarai^^, Issei Yamamura^, and Chisato Yamauchi^ 

^ Department of Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi, 464-860, Japan 
^ Department of Astronomy, Faculty of Science, University of Tokyo, 3-1-1 Hongo, Bunkyo-ku, Tokyo, 113-0003, Japan 
^ Institute of Space and Astronautical Science (ISAS), Japan Aerospace Exploration Agency (JAXA), 3-1-1, 
Yoshinodai, Sagamihara, Kanagawa, 229-8510, Japan 

* European Space Astronomy Center (ESAC), Villanueva de la Canada, Apartado 78, 28691 Madrid, Sapin 
^ INAF, Istituto di Fisica dello Spazio Interplanetario, Via del Fosso del Cavaliere 100, 00133 Roma, Italy 

® Dipartimento di Fisica, Universita' Roma Tre, Via della Vasca Navale 100, 00146 Roma, Italy 
^ Radio Astronomy Laboratory, University of California, Berkeley, USA, 

* Department of Physics, Faculty of Science, University of Tokyo, 3-1-1 Hongo, Bunkyo-ku, Tokyo, 113-0003, Japan 
® National Astronomical Observatory of Japan, Mitaka, Tokyo, 181-8588, Japan, 

Max-Planck-Institut fiir extraterrestrische Physik, Giessenbachstra/3e, 85748 Garching, Germany 
Academia Sinica, Institute of Astronomy and Astrophysics (ASIAA), Taipei 10617, Taiwan 

Korea Astronomy and Space Science Institute (KASI), 61-1, Hwaam-dong, Yuseong-gu, Daejeon, 305-348, Republic 
of Korea 

Graduate School of Science, Osaka University, 1-1, Machikaneyama, Toyonaka, Osaka, 560-0043, Japan 
Institute of Astronomy, Faculty of Sicence, University of Tokyo, Mitaka, Tokyo, 181-8588, Japan 
Space Applications Mission Directorate, Japan Aerospace Exploration Agency (JAXA), 2-1-1, Sengen, Tsukuba, 
Ibaraki, 305-8505, Japan 

Preprint online version: March 1, 2010 

ABSTRACT 

Context. AKARI is the first Japanese astronomical satellite dedicated to infrared astronomy. One of the main purposes 
of AKARI is the all-sky survey performed with six infrared bands between 9 and 200 /im during the period from 2006 
May 6 to 2007 August 28. In this paper, we present the mid-infrared part (9 /im and 18 /im bands) of the survey carried 
out with one of the on-board instruments, the Infrared Camera (IRC). 

Aims. We present unprecedented observational results of the 9 and 18 nm AKARI all-sky survey and detail the operation 
and data processing leading to the point source detection and measurements. 

Methods. The raw data are processed to produce small images for every scan and point sources candidates, above the 
5a"noise level per single scan, are derived. The celestial coordinates and fluxes of the events are determined statistically 
and the reliability of their detections is secured through multiple detections of the same source within milli-seconds, 
hours, and months from each other. 

Results. The sky coverage is more than 90% for both bands. A total of 877,091 sources (851,189 for 9/im, 195,893 for 
18 /xm) are confirmed and included in the current release of the point source catalogue. The detection limit for point 
sources is 50mJy and 90mJy for the 9 fim and 18 /im bands, respectively. The position accuracy is estimated to be 
better than 2". Uncertainties in the in-flight absolute flux calibration are estimated to be 3% for the 9 fim band and 
4% for the 18 /im band. The coordinates and fluxes of detected sources in this survey are also compared with those of 
the IRAS survey and found to be statistically consistent. 

Key words, infrared: general, galaxy: general, surveys, methods: observational, techniques: image processing 
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1. Introduction 

Unbiased and sensitive all-sky surveys at infrared wave- 
lengths are important for the various fields of astronomy. 
The first extensive survey of the mid- to far- infrared sky 
was made by the Infrared Astronomy Satellite (IRAS) mis- 
sion launched in 1983 (Neugebauer et al., 1984). IRAS sur- 
veyed 87% of the sky in four photometric bands at 12, 
25, 60 and 100 fiui and substantially pioneered the vari- 
ous new fields of astronomy, such as circumstellar debris 
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disks around Vega-like stars (Aumann et al., 1984), and a 
new class of galaxies that radiate most of their energy in 
the infrared (Soifer et al., 1987). 

A decade later than IRAS, the Midcoursc Space 
Experiment (MSX; Price et al., 2001) surveyed the Galactic 
plane as well as the regions not observed by or confused in 
the IRAS mission with higher sensitivity and higher spatial 
resolution (18.3") in four mid-infrared broad bands cen- 
tered at 8.28, 12.13, 14.65 and 21.23 /im and two narrow 
bands at 4.29 and 4.35 /um. The MSX catalogue (version 
1.2) of the Galactic plane survey contains 323,052 sources, 
3 times as many as IRAS listed for the same region. 

In 2006, AKARI, the first Japanese space mission ded- 
icated to infrared astronomical observations (Murakami 
et al., 2007), was launched and brought into a sun- 
synchronous polar orbit at an altitude of 700 km. It has two 
scientific instruments, the Infrared Camera (IRC; Onaka et 
al., 2007) for 2 26 fim and the Far-Infrared Surveyor (FIS; 
Kawada et al., 2007) for 50-200 /xm. AKARIhas a Ritchey- 
Chretien type cooled telescope with a primary-mirror aper- 
ture size of 685 mm (Kaneda et al., 2007), which is operated 
at 6 K by liquid helium and mechanical coolers. One of the 
major observational objectives of AKARI is an all-sky sur- 
vey observation. The survey was executed during the life 
time of the cooling medium between 2006 May 8 and 2007 
August 28. The 9 and 18 /xm bands of the IRC and the 60, 
90, 140, and 160 yum bands of the FIS were used for the 
all-sky survey. 

In this paper, we present the mid- infrared part of the all- 
sky survey performed with the IRC. The IRC was originally 
designed for imaging and spectroscopic observations in the 
pointing mode, however, the all-sky observation mode was 
added as an operation mode following ground tests, in 
which the acceptable performance of continuous survey- 
type observations was confirmed (Ishihara et al., 2006a). 
The data of the IRC all-sky survey observation have been 
processed by a dedicated program and a point source cat- 
alog has been prepared. The content of this paper is based 
on the 13-1 version of the AKARI/IRC All-Sky Survey Point 
Source Catalogue. 

The outline of the observation is presented in §2. The 
data reduction is described in §3. The quality of the cat- 
alogue is statistically evaluated in §4, and a summary is 
given in §5. 



2. Observations 

2.1. The AKARI satellite 

AKARI has two observational modes: the all-sky survey 
and pointed observations. In the pointed observations, the 
telescope stares at the target or makes round trip scans 
around the target for about 10 minutes. In the all-sky sur- 
vey, the spacecraft spins around the Sun-pointed axis once 
every orbit keeping the telescope toward a great circle and 
making continuous scans of the sky at a scan rate of 216" 
s~^. The orbit rotates around the axis of the Earth at the 
rate of the yearly round of the earth. Thus, the whole sky 
is in principle covered in half a year. 

In the 1^* half year (Pliasc 1) of the mission the all- 
sky survey dedicated as a first priority with pointed ob- 
servations toward the North Ecliptic Pole and the Large 
Magellanic Cloud. In the 2"'^ and 3'''^ half year (Phase 2a 
and 2b), the all-sky survey was continued until exhaust 



Table 1. Parameters for Mid-Infrared All-Sky Survey operation. 
Filter band (Camera) S9W (MIR-S) L18W (MIR-L) 



Wavelength* 
Isophotal wavelength 
Effective bandwidtii 
Sampling rate (period) 
Scan rate (exposure*) 
Reset rate (period) 
Operation 
Operated row 
Binning 

Virtual pixel scale 
Detection limit {5a) 



6.7-11.6 /im 13.9-25.6 /xm 
8.61 urn 18.39 /xm 

4.10 /im 9.97 /xm 

22.27 Hz (44 ms) 
216"s^i (11ms) 
0.074 Hz (13.464 s) 
256 X 2 pix 
117"", 125*'' 
4x1 pix 
9'.'36 X 9'.'36 10"4 x 9'.'36 
50 mJyt 120 mjyt 



* Effective exposures for point sources are determined not by the 
sampling rate but by the dwelling time of a source on a pixel, 
t Estimated value from readout noise in shuttered configuration 

in-orbit. 

J Defined as where the responsivity for a given energy is larger 
than 1/e of the peak. 



of the cooling medium (liquid helium). The time was dev- 
ided between pointed observations survey observations to 
increase the final sky coverage of the all-sky survey. 

2.2. The Infrared Camera (IRC) 

The mid-infrared component of the AKARI All-Sky Survey 
was performed with one of the two focal-plane instruments: 
the Infrared Camera (IRC). IRC covers the wavelength 
range 2-26 /xm with three independent channels: NIR (2- 
5.5 /im), MIR-S (6 12 /im) and MIR-L (12-26 /im). IRC was 
primarily designed for deep imaging and spectroscopy in 
pointed observations. All the channels have filter wheels 
which hold 3 filters and 2 spectroscopic dispersers. Each 
channel has a large format array that provides a wide field- 
of-view (FOV) of 10' x 10'. The MIR-S and MIR-L channels 
have infrared sensor arrays of 256 x 256 pixels (Si:As/CRC- 
744 manufactured by Raytheon) . The pixel scales for MIR-S 
and MIR-L are 2'.'34 x 2'.'34 and 2"51 x 2"39, respectively. 
The field-of-vicws of MIR-S and MIR-L are separated by 
20' in the cross-scan direction. More details and in-flight 
performance are described in Onaka et al. (2007). 

2.3. All-Sky Survey operations of the IRC 

Array Operation During the all-sky survey observations, 
only two out of 256 rows in the sensor array are operated 
in the continuous and non-destructive readout mode (scan 
operation of the array; Ishihara et al., 2006a). The 1^* row 
used in the operation can be selected arbitrarily out of 256 
rows and we adopted 117th row (hereafter row#l). The 2"*^ 
row (hereafter row:^2) was fixed to be 8 rows from the I*'* 
row the sampling rate of the array and the designed scan 
rate of the satellite. 

The sampling rate was set to 22.72 Hz (one sampling 
per 44 ms) taking account of the array operation conditions 
and the data down-link capacity. All the pixels are reset at 
a rate of 0.074 Hz (one reset per 306 samplings) to discharge 
the photo-current. 

The NIR channel is not used during the all-sky survey 
because of the capacity of the down link rate and because 
the alignment of the NIR array is not suited for the all-sky 
survey observation. 
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Fig. 1. Relative spectral response curve of the S9W (solid 
curve) and L18W bands (dashed) in units of electron/energy 
normalized to the peak. The system response curves 
of the IRAS 12 /im (dotted) and 25 /^m bands (dashed- 
dotted) are also shown for comparison. The RSRs of 
the IRC are available at http://www.ir.isas.jaxa.jp/ASTRO- 
F/Observation/RSRF/IRC_FAD/index.html. The system re- 
sponse curves of the IRAS are taken from Table II. C. 5 of IRAS 
Explanatory Supplement (Beichman et al., 1988). 



Exposure time The scan speed of the satellite in the survey 
observation mode is 216"s~^. As the pixel scale of the de- 
tector in the in-scan direction is about 2'.'34, the resulting 
effective exposure time for a point source is 11ms. 

Filter bands The all-sky survey is performed with two broad 
bands centered at 9 ^m (S9W filter of MIR-S) and 18 ^m 
(L18W filter of MIR-L). The relative spectral response 
(RSR) curves of the two bands are shown in Fig. 1 together 
with those of the IRAS 12 and 25 /im bands. 

Pixel scale The effective pixel size in the survey observation 
(hereafter virtual pixel size) is 9'/36 x 9'/36 for the 9/im 
band, and 10'.'4 x 9'.'36 for the 18 ^m band. The in-scan 
pixel size is fixed by the sampling rate, while the outputs of 
4 neighboring pixels are coadded to meet the down-link rate 
requirements and still make effective observations. A finer 
resolution is obtained in the data processing by combining 
images produced by the two rows. A detailed description 
of the reconstruction of the image and the resulting spatial 
resolution are given in Appendix A. 

Confirmation strategy It is difficult to distinguish events of 
real celestial objects from those due to cosmic ray hits only 
from the shape of signal because the virtual pixel size is not 
small enough to sample the PSF. We improve the reliability 
of source detections by adopting a three-step confirmation 
scheme. 

The detection of a celestial source in row#l is confirmed 
by a 2"'^ detection in row#2 87milli-seconds later. In this 
process, a large fraction of the signals due to cosmic ray hits 
are expected to be removed (milli-seconds confirmation). 

The width of the array in the cross-scan direction is 
about 10'. The scan path shifts at most by 4' (on the eclip- 
tic plane) due to the orbital motion. Thus, the scan path 
overlaps at least by 6' with the next scan 100 minutes after. 




Fig. 2. All-Sky image taken with the AKARI 9/xm band, from 
which the zodiacal light is subtracted simply as the low spatial 
frequency component derived from the raw data. 



The detection is thus comfirmed by the next scan osberva- 
tion. (hours confirmation). 

Furthermore, the scan path rotates by 180 degrees 
around the axis of the earth in half a year, giving another 
chance of detection six months later (months confirmation) . 
Objects, such as asteroids, comets and geostationary satel- 
lites, can be distinguished from stars and galaxies by hours- 
and/or months- confirmations. 

Sky coverage Large portions of the sky have a chance to 
be covered more than three times during the survey period 
(Phasel, Phase2a, and Phase2b; §2.1). During the semi- 
continuous survey, observational gaps (where no data are 
effectively available) appear because: (1) The survey ob- 
servations are halted for every pointed observation. A sky 
area of about 10' x 120° is skipped during a pointed obser- 
vation. (2) During Phase 2 the attitude of the satellite was 
operated actively also in the all-sky survey mode to cover 
the sky area not observed in Phase 1 in order to complete 
the far- infrared FIS all-sky survey (offset survey) . However 
it is the intermission of the continuous survey for the mid- 
infrared channels which have FOV directions different from 
the FIS. (3) The mid-infrared survey observations are some- 
times halted because of the limited downlink capacity. (4) 
The attitude of the satellite is lost due to a star tracking 
failure. (5) The shutter is closed in the moon avoidance re- 
gion. (6) The data taken during resets and during times 
affected by the heavy reset anomaly are masked by the 
pipeline software (§3.1). (7) Saturated pixels can not ob- 
serve the sky until having been reset (saturation trail). 

Taking into account (1), (2), and (3), which are auto- 
matically detected from the telemetry data, the final sky 
coverage is higher than 90% for both bands. Note that the 
actual sky area in which the source confirmation is carried 
out is smaller than this number (see §3.5). Fig. 2 shows the 
AKARI'd^ia low- resolution intensity map derived through 
the data processing described in this paper, from which the 
zodiacal light has been subtracted. 



3. Data processing 

The outline of the data processing (Fig. 3) for deriving 
the first point source catalogue of the AKARI/lYlC mid- 
infrared all-sky survey is summarized in the following: 
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— The raw telemetry data, which include the output of the 
sensor array, the house keeping data, and the output of 
the attitude control system, are down-linked from the 
satellite. Their timing is matched with each other and 
then they are registered to the database. 

— The raw data that are sandwiched by two successive 
resets in the scan are processed to make pieces of images 
corresponding to 10' x 50' wide sky regions (hereafter 
unit images) (Basic process; §3.1). 

— Then signals of the point source detections (events) are 
extracted from each processed image (Event detection 
and milli-seconds confirmation; §3.2). 

~ Next, the coordinates of all the events are determined 
according to the output of the attitude control sys- 
tem and refined using cross-correlations between the 
detected events and the prepared standard stars for the 
position determination (Pointing reconstruction; §3.3). 

— After the absolute flux calibration based on the mea- 
surements of standard stars, the fluxes of all the de- 
tected events are statistically derived (Flux calibration; 
§3.4). 

— Finally, reliable events are compiled into point source 
lists and the source lists of the 9 /xm band and 18 fim 
band are merged together to produce the IRC all-sky 
survey catalog (Catalogue compilation; §3.5). 

3.1. Basic process 

The basic process derives a unit image of 10' x 50' from the 
raw data for the extraction of point source candidates. 

ADU to electron conversion First, the pixel value of the raw 
data is converted from ADU into electrons as 



1.9e+05 



(1) 



where i is the pixel number, t is the time from the lat- 
est reset, S'("^(z,t) is the pixel value at {i,t) as the result 
of applying the nth step, and CF is the conversion factor 
measured from laboratory tests and assumed to be a con- 
stant for the pixels in the detector array (Ishihara et al., 
2003). 

Reset anomaly correction An anomalous behavior of the 
output level that persists a few seconds after the reset (re- 
set anomaly) is corrected. An example of this phenomenon 
and its correction are shown in Fig. 4. The offset level of 
the output of the detector is fairly sensitive to the tem- 
perature (Ishihara et al., 2003). We thus suppose that this 
phenomenon is explained by the drift of the offset level of 
the read-out circuit hybridized to the detector array owing 
to the temperature drift invoked by the reset current to 
discharge the stacked photo-electrons. Assuming that this 
behavior is pixel-independent in the detector array, the re- 
set anomaly is corrected as 



(2) 



where i is the pixel number, t is the time from latest reset, 
Cres is (pixel average of) the amount of the photo-electrons 
discharged during the latest reset, and 0{t, e^es) is the offset 
level variation. The offset level variation 0{t,eres) due to 
the reset is represented in the form of 
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Fig. 4. Example of the reset anomaly correction. The output 
signal level of a pixel in units of electron is plotted against the 
time from latest reset. This data set was taken under the con- 
stant illuminating source in the laboratory test. The raw output 
of the sensor (solid points) is fitted by the reset anomaly func- 
tion 0(t, Cint) (dotted curve) and corrected (dashed line). 

where Ci and C2 arc constants statistically derived from 
laboratory tests. 



Masking invalid data A masking pattern is created to ignore 
unusable data. (1) All the pixels are masked during the re- 
sets and the periods heavily affected by the reset anomaly. 
(2) Saturated pixels are masked from the moment of the 
saturation to the next reset. (3) The pixels under the slit 
masks are always masked. Both MIR-S and MIR-L channels 
have a slit and slit masks at the edge of the FOV for spec- 
troscopic observations of extended objects. The slit mask 
covers a few pixels located near the edge of the rows used 
in the survey operation. 

Linearity correction We correct for the non-linearity of the 
photo-response (L(eint)), which is supposedly due to the 
decrease in the bias voltage imposed on the detector array 
by the accumulation of the photo-electrons in each pixel. 
It is supposed to be a function of the number of stacked 
electrons (cint). An example of the non-linearity is shown 
in Fig. 5. The non-linearity of the photo-response is thus 
corrected as 



^(3)(^,t) = ^(2)(,^i)/i(ej„,), 



(4) 



assuming that their behaviors are the same for all the pix- 
els. The non-linearity function L(eint) is approximated by 
a spline function: 
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■^j(eint) = > ' C- 'e 



i=0 



(5) 



where j (j=0,l,2) represents the j range of Cint and each 



coefficient (C^*'') is derived from laboratory tests. 



Differentiation The signal is differentiated with respect to 
time. 



(6) 



0{t,eint) = Ciexp{~t/C2), 



(3) 



where At is the time step (44 ms) corresponding to the 
sampling rate. 
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Fig. 3. Outline of AKARI /IRC all-sky survey data processing. 
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Fig. 5. Example of the non-linearity of the photo-response of 
a pixel. Left panel shows time from reset versus output level 
of two pixels at the constantly illuminated condition in a lab- 
oratory test. The output profiles are fitted by linear functions 
using the data in the range, where the photo-response is still 
linear because the number of stacked electrons is small (7000- 
15000 ADU). Right panel shows the number of charged electrons 
(Y-axis) versus the deviation of the output level from the linear 
functions (X-axis). 



Flat fielding The differentiated signal S''^\i,t) is corrected 
for flat-fielding. 



(7) 



where F{i) is the normalized fiat correction factor for each 
pixel operated in the survey mode representing the disper- 
sion of the photo-response among the pixels. The flat func- 
tion F{i) is derived from multiple detections of the stars on 



different pixels. The background sky data are not used to 
avoid the effects of the scattered light. Fig. 6 shows the flat 
functions for both rows in both bands. 



Image construction from data of the two rows The data from 
row^l and row^2 are combined into a single image with 
a pixel scale of 1'.'56 to reject cosmic ray hits and produce 
a finer and higher S/N image. Details on the combining 
process are given in Appendix A. The mask patterns from 
the two rows are also combined in the same manner. Fig. 7 
shows examples of the raw data obtained by the two rows, 
the processed data of the two rows, (i.e. the resulting com- 
bined image of 10' x 50'), and the combined mask pattern. 
The images and mask patterns are stacked as the basic 
calibrated data and put into the next step, the source ex- 
traction process. 

3.2. Event detection and milli-seconds confirmation 

After the basic processing, signals above 5 a per scan are 
extracted from each image and checked by milli-seconds 
confirmation for point source detection (hereafter event). 
An example of these processes is shown in Fig. 8. 

First, events are extracted from two images obtained 
independently by the two rows, using the Source Extractor 
(Bertin et al., 2000) with a 3a threshold. Then events 
above Scr are extracted again on the finner combined image. 
Events from row#l, row#2 and the combined image are 
cross-identified by the detected position. Celestial sources 
are expected to be detected in both rows on the same pixel 
number (i) with similar fluxes, whereas false detections such 
as cosmic ray hits are expected to be detected only on one 
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Fig. 6. Flat functions for the one dimensional array for the band (left) and 18 fim band (right). The open circles represent 

row^^l and open squares represent row#2 in both panels. 




Fig. 7. Example of the processed data, (left) Raw data from 
Row#l and Row#2 are shown separately, (middle) Processed 
data. The data from both rows are combined, (right) Mask for 
the processed data. 

side. A pair of detections at the same sky position with 
similar fluxes (0.1 < Frow#i/^row#2 < 10) are selected as 
milli-seconds confirmed sources. Finally, the flux (photo- 
metric result) and the position (timing and pixel number 
of the detection) of the combined image are recorded for 
the confirmed events. 

3.3. Pointing reconstruction 

The celestial coordinate of a source (Qpixih t)) detected on 
the i-th pixel at the moment t are derived from 

Qobi{i,t) ^ Lit)Eii)QBit), (8) 

where Qsit) is the boresight position of the telescope at 
the moment t, E{i) is the conversion from the boresight to 
the ith pixel, and L{t) is the correction for the aberration 
due to the yearly revolution of the Earth and the orbital 
motion of the satellite. 

The improvement of boresignt position of the telescope 
(Qsit)) is carried out in collaboration with ESA by associ- 
ating detected events in the 9 and 18 /im bands and signals 
of the focal-star sensors at near-infrared on the focal plane 
with the stars in the positional reference catalogue prepared 




01 2345678 

distance (arcsec) 



Fig. 9. Statistical error of the pointing reconstruction using 50% 
of the catalogue data. The error is defined as the distance be- 
tween the position determined by the pointing reconstruction 
and the position of the input reference catalogue for events in 
the remaining 50% data. The lines show the fraction of the IRC 
events with an error smaller than given values. The color of the 
lines denotes the flux range of events. 



from MSX, 2MASS and IRAS (pointing reconstruction). 
The actual alignment of the FOV of each pixel E{i) includ- 
ing the effect of the distortion is optimized statistically in 
the pointing reconstruction process. 

Fig. 9 shows the accuracy of the pointing reconstruc- 
tion. The plot shows the fraction of the IRC events with 
an error smaller than given values. The error is estimated 
from the distance between the positions determined from 
the pointing reconstruction and those from the position 
reference catalogue. To make fair evaluation, the pointing 
reconstruction for this test is carried out using randomly 
selected sources amounting to half of the catalogue, and 
then the positions of the sources from the other half of the 
catalogue are determined for the evaluation. The error in- 
cludes the pointing reconstruction processing and the mea- 
surement errors. For the brightest sources the measurement 
error should be a minor contributor. We conclude that the 
position accuracy is better than 3" for 95 % of the events. 

Details of the pointing reconstruction are summarized in 
Salama et al. (2010). We adopt the position of each detected 
event derived from the pointing reconstruction process. 
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3.4. Flux calibration 

The photometric output for all the milli-seconds confirmed 
events is converted from electrons into physical units (Jy). 
The conversion function from electrons to Jy is derived sta- 
tistically by comparing the model fluxes with the measure- 
ments of hundreds of standard stars. The standard stars are 
selected from the infrared standard star network consisting 
of K- and M-giants (Cohen et al., 1999) and additional faint 
standard stars located around the north and south ecliptic 
poles (Reach et al., 2005; Ishihara et al., 2006b), which 
have a high visibility for the AKARI survey. The expected 
fluxes (/^"°*'^'^(Ai)) of the standard stars at the effective 
wavelengths (9 and 18 /xm) are calculated for the incident 
spectrum of f\ oc by convolving the model spectra 
of the standard stars (/a (A)) with the relative spectral re- 
sponse curves (Ri) in electron units as 



/'quoted / \ \ 



J R,iX)Xfxi\)dX 
J{^)R,{X)Xd\ ' 



(9) 



where i represents the band i, Xi is the effective wavelength 
of the band. 

The fitting function to convert the measured signals P 
into fluxes F is given by 
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Fig. 11. Ratio of the measured fluxes to the average fluxes as a 
function of observing time for five bright (> 1 Jy) standard stars 
observed at 9/im more than 30 times during the survey. The 
stars used and the corresponding 9 fluxes and number of de- 
tections are: HD42540 (plus, 8.33 Jy, 44 times), HD45669 (cross, 
9.79 Jy, 43 times), HD53501 (star, 9.30 Jy, 45 times), HD55865 
(open box, 16.7 Jy, 33 times) and HD170693 (filled box, 9.60 Jy, 
76 times), respectively. The closeup shown in the right bottom 
panel is to check for the temporal variation on a short timescale. 



log(F) = ^aiog(P)\ 



(10) 



i=0 



where C^'s are fitting coefficients (their values are given in 
the release note). Finally, the derived conversion functions 
are applied to all the milli-seconds confirmed events. The 
zero magnitude flux is 54.5 ± 0.845 Jy and 12.3 ± 0.204 Jy 
for the 9 /im and 18 /J,m bands, respectively. 

Fig. 10 (top) shows the flux derived (by Eq. 10) from 
the pipeline output as a function of the predicted in-band 
flux of the standard stars. 

The accuracy of the calibrated flux is investigated by 
the ratio of the measured flux to the predicted flux of the 
standard stars. The results shown in the bottom panels of 
Fig. 10 indicate an accuracy for the absolute calibration of 
about 3% for the 9 /xm band and 4% for the 18 ^m band. 

A systematic offset between A-type and K-M type giant 
standard stars is reported by Reach et al. (2005) for the 
calibration of the IRAC on Spitzer. We have only three A- 
type standard stars in our calibration for the 9 /xm band and 
no A-type stars are used for the calibration of the 18 fim 
band due to the detection limit. No systematic offset is 
seen in the present calibration between A-type stars and K- 
M giants within the measurement uncertainties (Fig. 10). 
Note that the calibration of pointing observations of the 
IRC imaging mode does not show any appreciable offsets 
either (Tanabe et al., 2008). 

To test of the long-term stability of the photo response 
we investigate flve standard stars that have been observed 
more than 30 times during the time of the survey. Fig. 11 
shows the ratio of the fluxes of individual measurements to 
the average fluxes of these stars as a funcion of time. From 
these data we deduce that the sensitivity is stable at the 
~2% level during the entire period of the observations. 

Laboratory measurements of the filter transmission in- 
dicate possible blue leaks between 3 and 4 /im for the 9 /im 
band and between 6 and 7^m for the 18 /im band. They 
are 0.01% at maximum and much smaller than the mea- 
surement errors. Thus the presence of the blue leaks is not 



confirmed. We calculate the predicted fluxes of the stan- 
dard stars with and without the blue leak and confirmed 
that the difference is less than 0.1%. We have verified that 
the blue leak of the 18 /xm band is negligible (< 1%) com- 
pared to the systematic errors by using asteroid calibrators 
whose flux can be well predicted (Miiller & Hasegawa, pri- 
vate communication). 

3.5. Catalogue compilation 

The flnal list of the sources (hereafter Point Source 
Catalogue; PSC) has been prepared based on the follow- 
ing criteria. First, groups of multiple events located within 
a region of 5" radius are recognized as a same source. In 
this process, events in the South Atlantic Anomaly (SAA) 
are excluded. Only groups containing at least two events 
are recognized as an actual celestial source. After the first 
grouping process, the distance from a source candidate to 
the nearest one is investigated and if there are two or more 
groups within 7", we consider them as a single source asso- 
ciated with outskirts events. In this case we take the group 
with the maximum number of events as a source candi- 
date and discard the other groups. Then, the source lists 
in the 9 fim and 18 fim band are merged into a single list. 
Sources within 7" are regarded as the same source in both 
bands. The position (RA, DEC) and associated position 
error (the major and minor axes and the position angle) 
are calculated from the events in the 9 /im band only, if 
the number of available events is larger than or equal to 
2. Otherwise, these data are calculated from the events in 
the 18 /im band only. The flux and associated error of the 
source are estimated from the mean and the mean error 
of multiple measurements of the events, respectively. Event 
data near the edge of the image strips are excluded from 
the flux calculation unless the exclusion leaves only zero or 
one event. 

The numbers of the events obtained in each step of 
the process are summarized in Table 2. The spatial dis- 
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Table 2. Number of events at each step of the process for the 
preparation of the PSC. 





9 ^im 


18 fim 


Detected 


16,752,471 


9,655,059 


Milli-seconds confirmed 


4,929,586 


1,300,945 


Hours or Months confirmed 


851,189* 


195,893* 


Band merged 


877,091* 



* - Current values. The number of sources may be updated in 
further catalogue releases. 





the largest errors are associated with sources with the low- 
est fluxes. We find that the most probable error is 2-3 % 
and the probability of larger errors is small. The relative er- 
rors are smaller than 15% for 80% of the sources and smaller 
than 30% for 96 % of the sources. It should be stressed that 
at this stage we cannot quantify how much variable objects 
"artificially" contribute to an increase of the scatter in the 
observed fluxes. The inclusion of mid-IR variability of the 
objects wiU be considered in future catalogue releases. 

Fig. 15 shows a typical signal-to-noise ratio (S/N) as- 
sociated with the different flux levels. The figure indicates 
that the expected S/N is ~ 6 for the faintest sources of ^ 
0.045 Jy in the 9/im band, while it is ^ 3 for 0.06 Jy 
sources in the 18 fj,m band. The figure also shows that the 
S/N increases with increasing fiux, but above 0.6 Jy and 
0.9 Jy it becomes constant or slightly decreases to ^ 20 
and ~ 15 for 9 /im band and 18 /xm band, respectively. This 
leveling-off of the S/N is due to uncertainties in the data 
reduction process. Errors in the correction for the detector 
reset anomaly, linearity, and flat-fielding limit accuracy of 
the of the flux measurement. 



Fig. 12. Spatial distribution in Aitoff Galactic coordinates pro- 
jection of the 18 fj,m events rejected in the catalogue compilation 
process. The grids indicate the equatorial coordinates. Most of 
the events correspond to high energy particles due to the South 
Atlantic Anomaly (SAA), or to moving objects such as aster- 
oids and comets (aligned around |/?| < 10°), to geostationary 
satellites (|Dec| 0°), and to retired geostationary satellites 
(|Dec| < 10°). The colour version will be available in the online 
version, in which the SAA events are plotted in cyan, and the 
remaining objects in red. 



tribution of the rejected events (those without hours and 
months confirmation) is shown in Fig. 12. Most of the re- 
jected events are ascribed to asteroids, geostationary satel- 
lites, and the high-energy particle hits in the SAA. These 
objects are recognized as such thanks to their spatial distri- 
bution and characteristics features. It should be noted that 
the actual sky coverage for this catalogue is smaller than 
the value quoted in §2.3; this is due to the severe condition 
of iVcvcnts > 2 as well as to the exclusion of all data affected 
by the SAA. 
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4. Evaluation of the catalogue 

4.1. Spatial distribution 

Fig. 13 shows the spatial distribution of the number density 
of the catalogued sources in the 9 and 18 /.tm bands after 
removal of the rejected events. One can easily appreciate 
the presence of highly populated regions, which correspond 
to the Galactic plane, the Large Magellanic Cloud (LMC), 
the Small Magellanic Cloud (SMC), as well as the nearby 
star forming regions such as p Oph, Orion and Taurus. Note 
that the 18 pm sources are strongly concentrated near the 
thick Galactic plane, whereas the 9 pm source counts fall 
off smoothly with the galactic latitude. 

4.2. Flux accuracy 

The histogram in Fig. 14 shows the relative flux errors in 
the 9 and 18 pm bands as a function of flux. As expected, 
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Fig. 15. Signal to noise ratio as a function of the 9 ^m (top) 
and 18 pm band flux (bottom) . 
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4.3. Position accuracy 

A test of the AKARI position accuracy is indicated in 
Fig. 16, which shows the angular separation between 
AKARIsouices and the nearest 2MASS source (Cutri et al., 
2003). The figure indicates that about 73% of the sources 
have an angular separation < 1", nearly 95% of the sources 
have a separation < 2". On average, the mean separa- 
tion between AKARI and 2MASS coordinates of matching 
sources is 0.8±0.6". 
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Fig. 18. Histograms of the spatial extension of AKARI sources 
at 9 pirn and 18 /xm, (solid and dashed lines, respectively). The 
spatial extension is calculated as (a + b)/2, where a and b are 
the length in arcsecond for the major axis and the minor axes 
of the source image, respectivel. 



Angular Separation AKARI-MIR vs. 2MASS [arcsec] 



Fig. 16. Histogram of the angular separation between AKARI 
coordinates and the 2MASS coordinates for the common 
sources. 



4.4. Number of detections per source 

The histogram in Fig. 17 shows the number of detections 
per source, both in differential and integrated counts. As 
described in §3.5, all the catalogued sources are detected 
at least twice. Over 80% of the sources have more than 3 
detections and about half of the sources have more than 5 
detections (right side scale of Fig. 17). The average number 
of detections of the catalogued sources is 6. 



9um 




Ndet 

Fig. 17. Number of detections per source at 9 fim (solid lines) 
and 18 fim source (dotted lines) . The differential counts are 
shown as histograms in the scale on the left vertical axis and 
the integrated counts are shown by the lines with symbols in 
the scale on the right vertical axis. 



4.5. False detections and missing sources 

False detections False detections are in principle rejected in 
the two-step source confirmation process. However, there 
is an off-chance of including (1) high energy particle hit 
events, (2) slow moving objects, (3) ghosts produced in the 
camera optics, and (4) those triggered by the signal satu- 
ration. 



Missing bright sources Some of the bright sources which 
have IRAS measurements are not included in the cata- 
logue because they are (1) not recognized as a point source 
with the beam size of AKARI because of the spatial ex- 
tension (Fig. 18 shows histograms of the spatial size of 
the AKARI/MIR sources), (2) located in the area not cov- 
ered by this survey, or (3) observed only once. A total 1722 
IRAS sources with the good quality flag (/quaii2 = 2 and 
/quai25 = 3) are not included in both the 9/Ltm and 18 /zm 
sources. 



4.6. Source counts versus flux 

It is instructive to compare the overlall distribution of the 
source counts between the AKARI and IRAS PSCs. The 
histograms of the source counts in the two surveys are com- 
pared in Fig. 19. The large difference in terms of the lower 
cutoff flux in the source counts between the two surveys is 
a clear signature of the considerably higher sensitivity, up 
to a factor of nearly 10, of the AKARI survey compared to 
IRAS. 

The distribution of source counts of the subsets of the 
9 ^m sources with and without a 18 /xm counterpart are 
shown in the left panel of Fig. 20. Similarly that for the 
18 /im sources with and without a 9 /im counterpart is 
shown in the right panel of Fig. 20. Most of the 18 /im 
sources have 9 fiui detections, whereas fainter 9 /im sources 
do not have their 18 /im counter parts. 
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4.7. Completeness 

The completeness of a survey above a given flux level is 
usually defined as the fraction of true sources that can be 
detected above that level. It is difficult to apply this con- 
cept to the ^i(r.4i?/ survey because one should dispose of a 
statistically significant sample of true sources with known 
IRC fluxes. The standard stars used for the AKARI/IRG 
PSC calibration might not represent a statistically signifi- 
cant sample in view of the rather poor coverage at low flux 
levels. 

To assess the completeness of the AKARI-MXR survey 
we have thus taken a different approach based on the dis- 
tribution of sources according to their flux. Fig. 19 shows 
an interesting feature: source counts decline exponentially 
with increasing flux after a peak value is reached around 
0.1 Jy (9/im band) and 0.2 Jy (18 /xm band). One can thus 
make a reasonable assumption that the exponential decay is 
an intrinsic property of source counts at the relevant wave- 
lengths and, on this basis, one can define completeness as 
the ratio of the number of sources actually observed by the 
number of sources predicted by the above equations. The 
results of completeness are shown in Figure 21. From this 
figure we deduce the completeness ratios reported in Table 
3. 

To evaluate the detection limit of the survey, we can 
make use of the S/N characteristics shown in Fig. 15. ^From 
this figure one can deduce that, for S/N ~ 5, the detection 
limit is about 0.05 and 0.09 Jy in the 9 /im and 18 ^m bands, 
respectively. 

A summary of the completeness of the survey at various 
flux levels, together with the corresponding values of the 
S/N ratio are given in Table 3. These results are in fairly 
good agreement with a pre-launch prediction in Ishihara et 
al. (2006a). 

The completeness is also tested separately for different 
Galactic latitude ranges, |5| < 2.2°, 2.2° < |6| < 8.3°, and 
|6| > 8.3°, where the division is made to have roughly sim- 
ilar numbers of the objects. Fig. 21 indicates that the com- 
pleteness becomes worse in the Galactic plane {\h\ < 2.2 
degree), which can be attributed to either the source con- 
fusion or the spatially variable background. 



Table 3. Completeness and signal-to-noise ratio 



9 pim band 


Completeness 


5% 50% 


80% 


100% 


Flux[Jy] 


0.07 0.10 


0.11 


0.15 


S/N 


6.7 7.6 


8.5 


9.8 


18 /im band 


Completeness 


5% 50% 


80% 


100% 


Flux[Jy] 


0.10 0.16 


0.21 


0.28 


S/N 


5.0 6.5 


7.6 


8.7 



4.8. Comparison with the IRAS catalogue 

To further evaluate the reliability of the AKARI Rux mea- 
surements at 9/im and 18 /xm, we have carried out a com- 
parison of these with the IRAS 12 /im and 25 /im fluxes for 
common sources. 
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Fig. 21. Completeness ratio of the AKARI/IRC survey in the 
9 /tm and 18 fim bands. Plots are made also for each galactic 
regions, |6| < 2.2°, 2.2° < |6| < 8.3°, and |6| < 8.3°, separately. 



In Fig. 22 a comparison is given between AKARI and 
IRAS fluxes observations for the cross-identifled objects. It 
clearly appears that there is a close correlation between the 
two sets of data. This is especially true for stellar sources; 
indeed, many MIR sources such as compact H ii regions, 
reflection nebulae and planetary nebulae, have extended 
emission or are located within extended objects. In such 
cases, IRAS measures the total flux of the extended emis- 
sion because of the larger apertures (0f75 x 4f5 — 4.'6 pixel 
size), whereas AKARI measures the flux of the peak emis- 
sion on the extended objects with a smaller aperture (~ 9" 
beam size) as demonstrated in the case of a reflection neb- 
ula (Ishihara et al., 2007). Some sources are brighter in the 
IRAS PSC than in the AKARI PSC. The difference could 
be attributed to the effect of the difference in the spatial 
resolution. 



5. Summary 

The AKARI Mid-Infrared all-sky survey was performed 
with two mid- infrared broad bands centered at 9 and 18 /im. 
More than 90% of the entire sky was observed in both 
bands. A total of 877,091 sources (851,189 for 9 /im, 195,893 
for 18 /im) are detected and included in the present release 
of the point source catalogue. This AKARI Mid-Infrared 
point source catalogue is scheduled for the public release 
in 2010 after the prioritized period for the team members 
(Release Note, 2009). We present the spatial distribution, 
flux distribution, flux accuracy, position accuracy, and com- 
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pleteness of the sources in the AKARIMIR All-Sky Survey 
Catalogue version /3-1. 

The AKARI mid-infrared survey provides a unique 

data-set relevant to interstellar, circumstellar-, and extra 
galactic astronomy. This new deep, large survey is well 
suited to research in the various fields of astronomy, such as 
search for warm debris disks, of asteroid analogue (Fujiwara 
et al. 2009a, 2009b), and provides valuable information for 
the study of planet formation and other fields. The chemical 
compositions of an unprecedented number of dust forming 
asymptotic giant branch (AGB) stars can also be investi- 
gated (Ita et al., 2009) making use of the characteristics 
of the filter bands of this survey. Such studies enable new 
discussions on the structure, the star formation history and 
the cycle of matter in our Galaxy. The mid-infrared all-sky 
survey also provides a chance to detect highly obscured ac- 
tive galactic nuclei (AGNs). which are difficult to observe 
in previous optical or X-ray surveys. The population of ob- 
scured AGNs is expected to contribute to the hard X-ray 
background, which is not competely resolved into individ- 
ual sources yet Ueda et al. (2003). 
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Appendix A: Nest 4x4 mode operation and unit 
image construction 

Here we expand on the detail of the data acquisition in 
the all-sky survey mode and subsequent image reconstruc- 
tion method to describe the pixel scale (Fig. A.l). Though 
the original pixel scale of the MIR-S camera is 2'.'34 x 2'.'34 
(2'.'51 X 2'.'39 for MIR-L), we degrade the spatial resolution 
in the survey observation to reduce the amount of output 
data to meet the downlink capacity. The pixel scale of the 
data taken by a single row in the in-scan direction is ad- 
justed as 9'.'36 by setting the sampling rate as 22.7 Hz under 
the scan rate of the satellite of 216" s^^. The pixel scale in 
the cross-scan direction is adjusted as 9'.'36 by binning of 
the output of four adjacent pixels. The resulting effective 
pixel scale obtained in the scan by a single row (hereafter 
virtual pixel scale) is 4 times larger than the original pixel 
scale. The scan observation was made by two rows to enable 
two independent observations in the milli-seconds interval. 
The sampling timing and the combination of the binning 
pixels are adjusted to construct two independent grids on 
the sky. Because the scan rate was changed from the value 
expected before launch due to the change in the altitude 
from the planned 750 km to 700 km, the shift of the grids 
in the in-scan direction is not exactly a half of the virtual 
pixel size. 

Fig. A. 2 compares an image of the same star constructed 
from a single-row observation in the survey mode, that pro- 
cessed from a two-row observation in the survey mode, and 
that obtained in the imaging mode in the pointed observa- 
tion with the original pixel scale. The FWHM of the PSF 
for the 9 /^m and 18 /xm bands is 5'.'5 and 5'.'7, respectively. 
The virtual pixel size is larger than the PSF size. However 
the process of the two-row observations in the survey mode 
reconstructs the image with the spatial resolution compati- 
ble with that obtained in the pointed observations with the 
original pixel size (center of Fig. A. 2). 

Thus the process of the two-row observation allows us 
to make (1) confirmation of the source detection in milli- 
seconds after, (2) reduction in the output data rate, and 
(3) higher spatial resolution than that in the single-row 
operation. Apparent asymmetry seen in the pointed data 
(Fig. A. 2 (a)) can be attributed to the effect of the tele- 
scope truss. Part of the asymmetry seen in the survey data 
may also be attributed to the undersampling operation. No 
variation of the PSFs is recognized in the 9 /xm and 18 /xm 
bands during the liquid Helium period in the pointing mode 
data. 
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Row#1 Row#2 Combined 



Fig. 8. Example of the event detection and milli-seconds confirmation. The left image is obtained by row^^l and the center image 
is obtained by row#2. The right image is constructed by combining data from two rows. The crosses mark events extracted on 
a single row image but rejected in milli-seconds confirmation. The circles show events extracted on both rows with similar fluxes 
and thus confirmed. 
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Fig. 10. (Top left) Calibrated fluxes obtained from the pipeline output through Eq. (10) are plotted as a function of the predicted 
in-band fluxes of the standard stars in the 9 /im band. The slid line indicates y = x for reference. (Bottom left) Ratio of the 
measured to predicted fluxes for the the standard stars. The right panels show the same plots for the 18 /^m band. The symbols 
and lines are the same as in the left panels. The blue points show three A-type stars used for the 9 pm band calibration. All the 
standard stars used for the 18 ^m band calibration are K-M giants and no A-type stars are used. 





Fig. 13. Spatial number density distribution of the detected sources in the Galactic coordinates of the Aitoff projection. (Left) 
9 jim sources of 851,189 and (right) 18 ^im sources of 195,893. 
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Fig. 14. Relative flux errors (%) for 9 (left) and 18 fim sources (right) as a function of flux plotted in the density maps. 
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Fig. 19. The distribution of source counts as a function of the 9 fim and 18 fim flux for AKARI is compared with that from the 
IRAS survey at 12 ^m and 25 fim. 
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Fig. 20. Flux distributions of the detected sources. (Left) Flux distribution for the 9 fim sources with 18 fim detections (dashed), 
without 18 /im detections (dotted) and total (solid). (Right) Flux distribution for the 18 fxm sources with 9 /^m detections (dashed), 
without 9/im detections (dotted) and total (solid). 
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Fig. 22. (Left) Comparison of the AKARI 9 fim versus IRAS 12 /im fluxes for the cross-identified sources. The red dots indicate 
sources with /quaii2 = 3 in the IRAS PSC labeled as a star in the SIMBAD database. The blue dots shows the other sources with 
/quaii2=3 in the IRAS PSC. The green dots indicate IRAS FSC sources. (Right) Comparison of the AKARI 18 /xm versus IRAS 
25 fim fluxes for the cross-identifled sources. The symbols are the same as in the left panel. 




Fig. A.l. (Left) Illustration of the nest 4x4 mode array operation in the survey mode. (Right) Image reconstruction of the data 
taken in the nest 4x4 mode operation. 




Fig. A. 2. 9/im image of the standard star (HD 42525): (a) the full-resolution image obtained in the pointed observation, (b) the 
images processed from the single-row data in the survey observation, and (c) the images reconstructed from the two-row data in 
the survey observation. 



